This code is Script 4b for Kitchel et al. TITLE manuscript.

COPY FROM THIS SCRIPT

(We expect this to be the main analyses in the text)

SESSION INFO TO DO

Quitting from lines 19-49 (year_dissimilarity_BC_balanced_models_4b.Rmd) 
library(data.table)
library(vegan)
library(sf)
library(concaveman) #polygon around points
library(betapart) #allows us to partition beta diversity
library(geosphere)
library(ggpubr) #stat_regline_equation
library(nlme)
library(mgcv) #to make gam
library(cowplot)
library(lme4)

#Pull Dissimilarity Means
distances_dissimilarities_allyears.r <- readRDS(here::here("output","dissimilarities", "distances_dissimilarities_allyears.r.rds"))

#make survey and survey unit factors
distances_dissimilarities_allyears.r[,survey:=factor(survey)][,survey_unit:=factor(survey_unit)]

#adjust years
distances_dissimilarities_allyears.r[,year_adj := year-min(year)+1]

#add new variable for year in sequence per region
distances_dissimilarities_allyears.r[,first_year := min(year),.(survey_unit)]
distances_dissimilarities_allyears.r[,last_year := max(year),.(survey_unit)]

#distances_dissimilarities_allyears.r[,year_in_seq := year-first_year+1]

distances_dissimilarities_allyears.r[,years_sampled := last_year-first_year+1]

###Palette for Plotting Palette for plotting all 37 survey units

survey_unit.list <- levels(distances_dissimilarities_allyears.r[,survey_unit])

palette_42 <- c(
  "#5A5156", #AI
  "#DF00DB", #BITS-1
  "#DB8EDA", #BITS-4
  "#F6222E", #CHL
  "#F8A19F", #DFO-NF
  "#16FF32", #DFO-QCS
  "#325A9B", #EBS
  "#3283FE", #EVHOE
  "#FEAF16", #FR-CGFS
  "#fccb6d", #GMEX-Fall
  "#1C8356", #GMEX-Summer
  "#C4451C", #GOA
  "#85660D", #GRL-DE
  "#B0009F", #GSL-N
  "#BF79B8", #GSL-S
  "#1CBE4F", #ICE-GFS
  "#782AB6", #IE-IGFS
  "#90AD1C", #MEDITS
  "#6B003A", #NAM
  "#A75B00", #NEUS-Fall
  "#E3B072", #NEUS-Spring
  "#02E8B6", #NIGFS-1
  "#97E7D5", #NIGFS-4
  "#B00068", #Nor-BTS-3
  "#00B9E3", #NS-IBTS-1
  "#95E2F4", #NS-IBTS-3
  "#B3CE73", #NZ-CHAT
  "#689500", #NZ-ECSI
  "#364d02",#NZ-SUBA
  "#AAF400", #NZ-WCSI
  "#AA0DFE", #PT-IBTS
  "#7f9eb8", #ROCKALL
  "#FA0087", #S-GEORG
  "#DEA0FD", #SCS-Summer
  "#FCEF88", #SEUS-fall
  "#A59405", #SEUS-spring
  "#FCE100", #SEUS-summer
  "#544563", #SWC-IBTS-1
  "#a37fc7", #SWC-IBTS-4
  "#C075A6", #WCANN
  "#BDCDFF", #ZAF-ATL
  "#003EFF"  #ZAF-IND
)

color_link <- data.table(survey_unit = survey_unit.list,hex = palette_42)

Add names for plotting


name_helper <- data.table(Survey_Name_Season = c("Aleutian Islands",
                                    "Baltic Sea Q1",
                                    "Baltic Sea Q4",
                                    "Chile",
                                    "Newfoundland",
                                    "Queen Charlotte Sound",
                                    "Eastern Bering Sea",
                                    "Bay of Biscay",
                                    "English Channel",
                                    "Gulf of Mexico Summer",
                                    "Gulf of Alaska",
                                    "Greenland",
                                    "N Gulf of St. Lawrence",
                                    "S Gulf of St. Lawrence",
                                    "Iceland",
                                    "Irish Sea",
                                    "Mediterranean",
                                    "Namibia",
                                    "NE US Fall",
                                    "NE US Spring",
                                    "N Ireland Q1",
                                    "N Ireland Q4",
                                    "Barents Sea Norway Q3",
                                    "N Sea Q1",
                                    "N Sea Q3",
                                    "Chatham Rise NZ",
                                    "E Coast S Island NZ",
                                    "W Coast S Island NZ",
                                    "Portugal",
                                    "S Georgia",
                                  "Scotian Shelf Summer",
                                  "SE US Fall",
                                  "SE US Spring",
                                  "SE US Summer",
                                  "W Coast US",
                                  "Atlantic Ocean ZA",
                                  "Indian Ocean ZA",
                                   "Rockall Plateau",
                                  "Scotland Shelf Sea Q1",
                                  "Scotland Shelf Sea Q4",
                                  "Falkland Islands",
                                  "Gulf of Mexico Fall",
                                  "Barents Sea Norway Q1",
                                  "Sub-Arctic NZ",
                                  "Scotian Shelf Spring"),
                          survey_unit = c(
                                  "AI",        
                                  "BITS-1",    
                                  "BITS-4",    
                                  "CHL",       
                                  "DFO-NF",    
                                  "DFO-QCS",   
                                  "EBS",       
                                  "EVHOE",     
                                  "FR-CGFS",   
                                  "GMEX-Summer",
                                  "GOA",       
                                  "GRL-DE",    
                                  "GSL-N",     
                                  "GSL-S",     
                                  "ICE-GFS",   
                                  "IE-IGFS",   
                                  "MEDITS",    
                                  "NAM",       
                                  "NEUS-Fall", 
                                  "NEUS-Spring",
                                  "NIGFS-1",   
                                  "NIGFS-4",   
                                  "Nor-BTS-3", 
                                  "NS-IBTS-1", 
                                  "NS-IBTS-3", 
                                  "NZ-CHAT",   
                                  "NZ-ECSI",   
                                  "NZ-WCSI",   
                                  "PT-IBTS",   
                                  "S-GEORG",   
                                  "SCS-SUMMER",
                                  "SEUS-fall", 
                                  "SEUS-spring",
                                  "SEUS-summer",
                                  "WCANN",     
                                  "ZAF-ATL",   
                                  "ZAF-IND",
                                  "ROCKALL",
                                  "SWC-IBTS-1",
                                  "SWC-IBTS-4",
                                  "FALK",
                                  "GMEX-Fall",
                                  "Nor-BTS-1",
                                  "NZ-SUBA",
                                  "SCS-SPRING"
                          ))


color_link <- color_link[name_helper, on = "survey_unit"]

##Make GAMs

Bray Curtis

bray_curtis_balanced_gam <- gam(bray_curtis_dissimilarity_balanced_mean ~ year + s(survey_unit, year, bs = "fs", m = 1),#random smooth
                            data = distances_dissimilarities_allyears.r)

##Make LMERS

Bray These all converged

#running with lme instead of lmer gave same results, but allowed for calculation of p-value
bray_curtis_balanced_lme <- lme(bray_curtis_dissimilarity_balanced_mean ~ year_adj, random = (~1 + year_adj|survey_unit),data = distances_dissimilarities_allyears.r)

#but also run with lmer for confint
bray_curtis_balanced_lmer <- lmer(bray_curtis_dissimilarity_balanced_mean ~ year_adj + (1 + year_adj|survey_unit),data = distances_dissimilarities_allyears.r)

summary(bray_curtis_balanced_lme)
Linear mixed-effects model fit by REML
  Data: distances_dissimilarities_allyears.r 

Random effects:
 Formula: ~1 + year_adj | survey_unit
 Structure: General positive-definite, Log-Cholesky parametrization
            StdDev     Corr  
(Intercept) 0.17066632 (Intr)
year_adj    0.00259622 -0.769
Residual    0.03491895       

Fixed effects:  bray_curtis_dissimilarity_balanced_mean ~ year_adj 
 Correlation: 
         (Intr)
year_adj -0.79 

Standardized Within-Group Residuals:
        Min          Q1         Med          Q3         Max 
-4.11521011 -0.53411661  0.03450067  0.57016282  4.45172058 

Number of Observations: 897
Number of Groups: 42 
anova(bray_curtis_balanced_lme)

bray_curtis_balanced_coefs <- data.table(coef(bray_curtis_balanced_lme))
bray_curtis_balanced_coefs[,survey_unit := rownames(coef(bray_curtis_balanced_lme))][,Year := round(year_adj,5)][,Intercept := round(`(Intercept)`,2)]
View(bray_curtis_balanced_coefs)

bray_curtis_balanced_coefs <- bray_curtis_balanced_coefs[color_link, on = "survey_unit"]

bray_curtis_balanced_coefs.exp <- bray_curtis_balanced_coefs[,.(Survey_Name_Season, Intercept, Year)]

#export this table
fwrite(bray_curtis_balanced_coefs.exp, file = here::here("output","bray_curtis_balanced_coefs.exp.csv"))

Get LMER model as predictions


# need to sort out year in seq versus overall year models
#new data for lmer
lmer_year <- seq(min(distances_dissimilarities_allyears.r[,year]), max(distances_dissimilarities_allyears.r[,year]), by = 0.1)

lmer_year_adj <- seq(min(distances_dissimilarities_allyears.r[,year_adj]), max(distances_dissimilarities_allyears.r[,year_adj]), by = 0.1)

#predict average lmer
lmer_bray_balanced_predictions <- data.table(year = lmer_year, year_adj = lmer_year_adj)

#confidence intervals
bray_curtis_balanced_lmer_confint <- confint(bray_curtis_balanced_lmer)
Computing profile confidence intervals ...
#populate data table of lmer predictions
lmer_bray_balanced_predictions[,bray_curtis_lmer_preds := fixef(bray_curtis_balanced_lmer)[[1]] + fixef(bray_curtis_balanced_lmer)[[2]] * year_adj][,bray_curtis_lmer_preds_lowerCI := bray_curtis_balanced_lmer_confint[5] + bray_curtis_balanced_lmer_confint[6] * year_adj][,bray_curtis_lmer_preds_upperCI := bray_curtis_balanced_lmer_confint[11] + bray_curtis_balanced_lmer_confint[12] * year_adj]

###Move forward with Bray Curtis Balanced (others to supplement)

Coefficients for LMER by survey_unit

#unique survey unit year combos
survey_unit_sampling_years <- unique(distances_dissimilarities_allyears.r[,.(survey_unit, year_adj, year, years_sampled)])

# see group coefficients
model_coefs_reduced <- data.table(transform(as.data.frame(ranef(bray_curtis_balanced_lmer)), lwr = condval - 1.96*condsd, upr = condval + 1.96*condsd))
#https://stackoverflow.com/questions/69805532/extract-the-confidence-intervals-of-lmer-random-effects-plotted-with-dotplotra


#ONLY SLOPES
model_coefs_reduced <- model_coefs_reduced[term == "year_adj",]

model_coefs_reduced[,survey_unit := grp][,year_adj := condval]

#merge with duration of survey
model_coefs_reduced_length <- model_coefs_reduced[survey_unit_sampling_years, on = "survey_unit"]


model_coefs_reduced_length[,years_sampled := as.numeric(years_sampled)][,Directional_Change := ifelse(year_adj > 0, "Differentiation","Homogenization")]

#does it cross zero?
model_coefs_reduced_length[,significant := ifelse(lwr >0 & upr>0,T,ifelse(lwr<0 & upr<0,T,F))]

#significant directional change
model_coefs_reduced_length[,Directional_Change_sig := ifelse(year_adj > 0 & significant == T, "Differentiation",ifelse(year_adj < 0 & significant == T, "Homogenization", "No trend in\ndissimilarity"))]


#min max distances_dissimilarities
min_bray_reduced <- min(distances_dissimilarities_allyears.r[,bray_curtis_dissimilarity_balanced_mean], na.rm = T)
max_bray_reduced <- max(distances_dissimilarities_allyears.r[,bray_curtis_dissimilarity_balanced_mean], na.rm = T)

model_coefs_reduced_length <- model_coefs_reduced_length[color_link, on = "survey_unit"]

#delete any NAs
model_coefs_reduced_length <- na.omit(model_coefs_reduced_length, cols = "significant")

#order table by coefficient
setorder(model_coefs_reduced_length, year_adj)

BC_balanced_model_coefs_reduced_length.unique <- unique(model_coefs_reduced_length[,.(condval,condsd, lwr, upr, survey_unit, year_adj, years_sampled, Directional_Change, hex, Survey_Name_Season, significant, Directional_Change_sig)]) 

#extract color hexes
#year adj coef order
color_year_adj_order <- BC_balanced_model_coefs_reduced_length.unique[,hex]

#alphabetical order
BC_balanced_model_coefs_reduced_length.unique.alpha <- setorder(BC_balanced_model_coefs_reduced_length.unique, Survey_Name_Season)

BC_balanced_model_coefs_reduced_length.unique.alpha[,trend_color := ifelse(Directional_Change_sig == "Homogenization", "#e7ac5b", ifelse(Directional_Change_sig == "Differentiation","#91c874","#cbbfde"))]

color_alpha_order <- BC_balanced_model_coefs_reduced_length.unique.alpha[,hex]
color_alpha_order_bytrend <- BC_balanced_model_coefs_reduced_length.unique.alpha[, trend_color]

saveRDS(BC_balanced_model_coefs_reduced_length.unique, here::here("output","region_stats","BC_balanced_model_coefs_reduced_length.unique.Rds"))

Bar Plot Coefficient LMER

BC_BALANCED_Dissimilarity_Coef_errorbar_reduced <- ggplot() +
    geom_errorbar(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, ymin = lwr, ymax = upr), fill = "grey", width = 0) + #add confidence intervals
  geom_point(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, size = years_sampled, fill = Directional_Change_sig, color = Directional_Change_sig), stat = 'identity', shape = 21) +
  scale_fill_manual(values = c("white","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_color_manual(values = c("black","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_size_binned(range = c(1,8), name = "Survey period length") +
  geom_hline(yintercept = 0) +
  scale_y_continuous(breaks = seq(-0.005, 0.0075, by = 0.0025), labels = c("-0.005","","0", "", "0.005",  "")) +
  xlab("Survey unit") +
  ylab("β-diversity trend") + #balanced BC dissimilarity trend
  coord_flip() +
  theme_classic() +
  theme(axis.text.y = element_text(colour = color_year_adj_order, face = "bold"), axis.title.y = element_blank(), axis.text.x = element_text(size = 15), axis.title.x = element_text(size = 15), legend.position = c(0.25,0.8), legend.direction = "vertical")
Warning: Ignoring unknown parameters: `fill`Warning: Ignoring unknown aesthetics: labelWarning: Ignoring unknown aesthetics: labelWarning: Vectorized input to `element_text()` is not officially supported.
ℹ Results may be unexpected or may change in future versions of ggplot2.
#pull legend for homogenization
directional_change_legend_plot <- BC_BALANCED_Dissimilarity_Coef_errorbar_reduced + 
  scale_fill_manual(values = c("white","black","grey"), name = "Dissimilarity trend") +
  scale_color_manual(values = c("black","black","grey"), name = "Dissimilarity trend") +
  scale_size_binned(range = c(1,8), name = "Survey period length", guide = "none") +
  theme(legend.position = "right", legend.background = element_rect(fill= "transparent")) +
  guides(colour = guide_legend(override.aes = list(size=6)))
Scale for fill is already present.
Adding another scale for fill, which will replace the existing scale.Scale for colour is already present.
Adding another scale for colour, which will replace the existing scale.Scale for size is already present.
Adding another scale for size, which will replace the existing scale.
BC_BALANCED_Dissimilarity_Coef_errorbar_reduced


#ALT grey scale
BC_BALANCED_Dissimilarity_Coef_errorbar_reduced_greyscale <- ggplot() +
    geom_errorbar(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, ymin = lwr, ymax = upr), fill = "grey", width = 0) + #add confidence intervals
  geom_point(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, size = years_sampled, fill = Directional_Change_sig, color = Directional_Change_sig), stat = 'identity', shape = 21) +
  scale_fill_manual(values = c("white","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_color_manual(values = c("black","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_size_binned(range = c(1,8), name = "Survey period length") +
  geom_hline(yintercept = 0) +
  scale_y_continuous(breaks = seq(-0.005, 0.0075, by = 0.0025), labels = c("-0.005","","0", "", "0.005",  "")) +
  xlab("Survey unit") +
  ylab("β-diversity trend") + #balanced BC dissimilarity trend
  coord_flip() +
  theme_classic() +
  theme(axis.text.y = element_text(face = "bold"), axis.title.y = element_blank(), axis.text.x = element_text(size = 15), axis.title.x = element_text(size = 15), legend.position = c(0.25,0.8), legend.direction = "vertical")
Warning: Ignoring unknown parameters: `fill`Warning: Ignoring unknown aesthetics: labelWarning: Ignoring unknown aesthetics: label

Alternatively, we color this plot by trend experienced

#"#73BA4D","#E0962C","#cbbfde"

BC_balanced_Dissimilarity_Coef_errorbar_reduced_colorbytrend <- ggplot() +
    geom_errorbar(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, ymin = lwr, ymax = upr), fill = "grey", width = 0) + #add confidence intervals
  geom_point(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, size = years_sampled, color = Directional_Change_sig), stat = 'identity') +
  scale_color_manual(values = c("#73BA4D","#E0962C","#cbbfde"), name = "Dissimilarity trend", guide="none") +
  scale_size_binned(range = c(1,8), name = "Survey period length\n(years)") +
  geom_hline(yintercept = 0) +
  scale_y_continuous(breaks = seq(-0.005, 0.0075, by = 0.0025), labels = c("-0.005","","0", "", "0.005",  "")) +
  xlab("Survey unit") +
  ylab("β-diversity trend") + #balanced BC dissimilarity trend
  coord_flip() +
  theme_classic() +
  theme(axis.text.y = element_text(face = "bold"), axis.title.y = element_blank(), axis.text.x = element_text(size = 15), axis.title.x = element_text(size = 15), legend.position = c(0.3,0.8), legend.direction = "vertical", legend.text = element_text(size = 15), legend.title = element_text(size = 16))
Warning: Ignoring unknown parameters: `fill`Warning: Ignoring unknown aesthetics: labelWarning: Ignoring unknown aesthetics: label
directional_change_legend_plot_colorbytrend <- BC_balanced_Dissimilarity_Coef_errorbar_reduced_colorbytrend + 
  theme(legend.position = "right", legend.background = element_rect(fill= "transparent"), 
         legend.text = element_text(size = 15), legend.title = element_text(size = 16)) +
  guides(colour = guide_legend(override.aes = list(size=6)), size = "none")

Wavy Line Plot for GAMs

Generate predicted values


#add colors and names to full dissimilarity data table
distances_dissimilarities_allyears.r <- distances_dissimilarities_allyears.r[color_link, on = "survey_unit"]

#generate new data to smooth lines (need year and season survey combinations)
year_survey_unit_expand.dt <- data.table(survey_unit = as.character(NULL), year = as.numeric(NULL), year_adj = as.numeric(NULL ))

for (i in 1:length(survey_unit.list)) {
  #generate year vectors
  survey_unit_years <- unique(distances_dissimilarities_allyears.r[survey_unit == survey_unit.list[i],.(survey_unit, year, year_adj)])
  
  years <- seq(min(survey_unit_years[,year]), max(survey_unit_years[,year]), by = 0.1)
  
  year_adjust <- seq(min(survey_unit_years[,year_adj]), max(survey_unit_years[,year_adj]), by = 0.1)
  
  year_survey_unit_expand.dt_addition <- data.table(survey_unit = survey_unit.list[i], year = years, year_adj = year_adjust)
  
  year_survey_unit_expand.dt <- rbind(year_survey_unit_expand.dt, year_survey_unit_expand.dt_addition)
}

#add colors and names to full year and survey unit combination table
year_survey_unit_expand.dt <- year_survey_unit_expand.dt[color_link, on = "survey_unit"]

Get model as predictions

#for plotting, get model as predictions
bray_curtis_balanced_gam_predictions <- predict(bray_curtis_balanced_gam, se.fit = TRUE, newdata = year_survey_unit_expand.dt)

#merge into table
year_survey_unit_expand.dt[,bray_glm_mod_fit := bray_curtis_balanced_gam_predictions$fit][,bray_glm_mod_fit_SE := bray_curtis_balanced_gam_predictions$se.fit]

Produce Plot of GAM and mean LMER line

points_wavylines_bray_balanced_year_reduced_gam_nolmer <- ggplot() +
 # geom_ribbon(data = lmer_bray_balanced_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.2) +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r,cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_balanced_mean,
                 color = Survey_Name_Season), alpha = 0.5, size = 1) +
    geom_line(data = na.omit(year_survey_unit_expand.dt,cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season)) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt,cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
 # geom_line(data = lmer_bray_balanced_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  color_alpha_order, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])),
       y = c(0.15,0.9)) +
  xlab("Year") +
ylab("Balanced BC dissimilarity") +
  theme(legend.position = "null")

points_wavylines_bray_balanced_year_reduced_gam_nolmer

ggsave(points_wavylines_bray_balanced_year_reduced_gam_nolmer, path = here::here("figures"), filename ="points_wavylines_bray_balanced_year_reduced_gam_nolmer.jpg", height = 5, width = 5, unit = "in")


#with lmer

points_wavylines_bray_balanced_year_reduced_gam <- ggplot() +
  geom_ribbon(data = lmer_bray_balanced_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.3) +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r, cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_balanced_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
  geom_line(data = lmer_bray_balanced_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  color_alpha_order, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])),
       y = c(0.15,0.9)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_balanced_year_reduced_gam

ggsave(points_wavylines_bray_balanced_year_reduced_gam, path = here::here("figures"), filename ="points_wavylines_bray_balanced_year_reduced_gam.jpg", height = 6, width = 6, unit = "in")


#ALT
#plot all, but same color scheme (grey)
points_wavylines_bray_balanced_year_reduced_gam_greyscale <- ggplot() +
  geom_ribbon(data = lmer_bray_balanced_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.3) +
  geom_point(data = distances_dissimilarities_allyears.r,
             aes(x = year,
                 y = bray_curtis_dissimilarity_balanced_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = year_survey_unit_expand.dt,
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = year_survey_unit_expand.dt, aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
  geom_line(data = lmer_bray_balanced_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  rep("black", times = length(unique(distances_dissimilarities_allyears.r$Survey_Name_Season))), name = "Survey Unit") +
  scale_fill_manual(values =  rep("black", times = length(unique(distances_dissimilarities_allyears.r$Survey_Name_Season))), guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])),
       y = c(0.15,0.9)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_balanced_year_reduced_gam_greyscale

ggsave(points_wavylines_bray_balanced_year_reduced_gam_greyscale, path = here::here("figures"), filename ="points_wavylines_bray_balanced_year_reduced_gam_greyscale.jpg", height = 6, width = 6, unit = "in")

Alternative, color by trend


points_wavylines_bray_balanced_year_reduced_gam_colorbytrend <- ggplot() +
  geom_ribbon(data = lmer_bray_balanced_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.3) +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r, cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_balanced_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
  geom_line(data = lmer_bray_balanced_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  color_alpha_order_bytrend, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order_bytrend, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year]))) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_balanced_year_reduced_gam_colorbytrend

ggsave(points_wavylines_bray_balanced_year_reduced_gam_colorbytrend, path = here::here("figures"), filename ="points_wavylines_bray_balanced_year_reduced_gam_colorbytrend.jpg", height = 6, width = 6, unit = "in")

#plot each independently for supplement
#all survey names = 
all_survey_names <- sort(unique(distances_dissimilarities_allyears.r$Survey_Name_Season))
#list of plots
points_wavylines_bray_balanced_year_reduced_gam_individual <- list()
for (i in 1:length(all_survey_names)) {
points_wavylines_bray_balanced_year_reduced_gam_individual[[i]] <- ggplot() +
  geom_point(data = distances_dissimilarities_allyears.r[Survey_Name_Season == all_survey_names[i]],
             aes(x = year,
                 y = bray_curtis_dissimilarity_balanced_mean), alpha = 0.4, color = "black") +
    geom_line(data = year_survey_unit_expand.dt[Survey_Name_Season == all_survey_names[i]],
             aes(x = year,
                 y = bray_glm_mod_fit), alpha = 0.6) +
  geom_ribbon(data = year_survey_unit_expand.dt[Survey_Name_Season == all_survey_names[i]], aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE), alpha=0.1) + #add standard error
  theme_classic() +
#  lims(x = c(min(distances_dissimilarities_allyears.r[Survey_Name_Season == all_survey_names[i],year]),max(distances_dissimilarities_allyears.r[Survey_Name_Season == all_survey_names[i],year])),
#       y = c(0.15,0.9)) +
  xlab("Year") +
ylab("beta-diversity") +
  facet_wrap(~Survey_Name_Season, ncol = 5) +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

print(points_wavylines_bray_balanced_year_reduced_gam_individual[[i]])

}

saveRDS(points_wavylines_bray_balanced_year_reduced_gam_individual, here::here("figures","points_wavylines_bray_balanced_year_reduced_gam_individual.Rds"))

#print to pdf
library(gridExtra)

ggsave(
   filename = here::here("figures","points_wavylines_bray_balanced_year_reduced_gam_individual.pdf"), 
   plot = marrangeGrob(points_wavylines_bray_balanced_year_reduced_gam_individual, nrow=1, ncol=1), 
   width = 8.5, height = 11
)
Warning: Removed 1 rows containing missing values (`geom_point()`).Warning: Removed 1 row containing missing values (`geom_line()`).Warning: no non-missing arguments to max; returning -InfWarning: Removed 1 rows containing missing values (`geom_point()`).Warning: Removed 1 row containing missing values (`geom_line()`).Warning: no non-missing arguments to max; returning -InfWarning: Removed 1 rows containing missing values (`geom_point()`).Warning: Removed 1 row containing missing values (`geom_line()`).Warning: no non-missing arguments to max; returning -Inf

Merge BC versus Year plot with GAMS and Region vs. coefficient plot for LMERs


BC_BALANCED_GAM_LMER_merge_legend <- ggdraw(xlim = c(0, 40.5), ylim = c(0, 21)) +
    draw_plot(points_wavylines_bray_balanced_year_reduced_gam,
                                         x = 1, y = 1, width = 20, height = 20) +
    draw_plot(BC_BALANCED_Dissimilarity_Coef_errorbar_reduced +
        theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=16), #change legend title font size
        legend.text = element_text(size=14)), #change legend text font size),
                                         x = 20, y = 1, width =19, height = 20) +
    draw_plot(get_legend(directional_change_legend_plot + 
      theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=15), #change legend title font size
        legend.text = element_text(size=13))), #change legend text font size)
                                       x = 27, y = 12, width = 3, height = 2) +
  geom_text(aes(x = 2, y = 20.7), label = ("a."), size =8, fontface = "bold") +
  geom_text(aes(x =20, y = 20.7), label = ("b."), size =8, fontface = "bold")


ggsave(BC_BALANCED_GAM_LMER_merge_legend, path = here::here("figures"), filename = "BC_BALANCED_GAM_LMER_merge_legend.png", height = 8, width = 14, units = "in")

#ALT GREY SCALE
BC_BALANCED_GAM_LMER_merge_legend_greyscale <- ggdraw(xlim = c(0, 40.5), ylim = c(0, 21)) +
    draw_plot(points_wavylines_bray_balanced_year_reduced_gam_greyscale,
                                         x = 1, y = 1, width = 20, height = 20) +
    draw_plot(BC_BALANCED_Dissimilarity_Coef_errorbar_reduced_greyscale +
        theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=16), #change legend title font size
        legend.text = element_text(size=14)), #change legend text font size),
                                         x = 20, y = 1, width = 19, height = 20) +
    draw_plot(get_legend(directional_change_legend_plot + 
      theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=15), #change legend title font size
        legend.text = element_text(size=13))), #change legend text font size)
                                x = 27, y = 12, width = 3, height = 2) +
  geom_text(aes(x = 2, y = 20.7), label = ("a."), size =8, fontface = "bold") +
  geom_text(aes(x =20, y = 20.7), label = ("b."), size =8, fontface = "bold")

ggsave(BC_BALANCED_GAM_LMER_merge_legend_greyscale, path = here::here("figures"), filename = "BC_BALANCED_GAM_LMER_merge_legend_greyscale.png", height = 8, width = 14, units = "in")

#ALT COLOR BY TREND
BC_balanced_GAM_LMER_merge_legend_colorbytrend <- ggdraw(xlim = c(0, 40.5), ylim = c(0, 21)) +
    draw_plot(points_wavylines_bray_balanced_year_reduced_gam_colorbytrend,
                                         x = 1, y = 1, width = 20, height = 20) +
    draw_plot(BC_balanced_Dissimilarity_Coef_errorbar_reduced_colorbytrend +
        theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
       # legend.title = element_text(size=16), #change legend title font size
       # legend.text = element_text(size=14)
       ), #change legend text font size),
                                         x = 20, y = 1, width = 19, height = 20) +
    draw_plot(get_legend(directional_change_legend_plot_colorbytrend + 
      theme(legend.key.size = unit(0.5, 'cm'), #change legend key size
        legend.title = element_text(size=16), #change legend title font size
        legend.text = element_text(size=15))), #change legend text font size)
                                x = 27, y = 8, width = 3, height = 2) +
  geom_text(aes(x = 2, y = 20.7), label = ("a."), size =8, fontface = "bold") +
  geom_text(aes(x =20, y = 20.7), label = ("b."), size =8, fontface = "bold")

ggsave(BC_balanced_GAM_LMER_merge_legend_colorbytrend, path = here::here("figures"), filename = "BC_balanced_GAM_LMER_merge_legend_colorbytrend.png", height = 8, width = 14, units = "in")
LS0tCnRpdGxlOiAiWWVhciBCYWxhbmNlZCBEaXNzaW1pbGFyaXR5IE1vZGVscyAocmVwbGFjZW1lbnQvdHVybm92ZXIpIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCgpUaGlzIGNvZGUgaXMgU2NyaXB0IDRiIGZvciBLaXRjaGVsIGV0IGFsLiBUSVRMRSBtYW51c2NyaXB0LgoKQ09QWSBGUk9NIFRISVMgU0NSSVBUCgooV2UgZXhwZWN0IHRoaXMgdG8gYmUgdGhlIG1haW4gYW5hbHlzZXMgaW4gdGhlIHRleHQpCgotIFRoaXMgcHJvamVjdCBpcyBhIGNvbGxhYm9yYXRpdmUgZWZmb3J0IHRvIGRlc2NyaWJlIGNoYW5nZXMgaW4gdGF4b25vbWljIGNvbXBvc2l0aW9uICBvZiBmaXNoIGNvbW11bml0aWVzIGFyb3VuZCB0aGUgd29ybGQtLWFzIHNhbXBsZWQgYnkgYm90dG9tIHRyYXdsIHN1cnZleXMuCgotIENvZGUgYnkgWm/DqyBKLiBLaXRjaGVsCgpTRVNTSU9OIElORk8gVE8gRE8KCmBgYHtyIHNldHVwfQpsaWJyYXJ5KGRhdGEudGFibGUpCmxpYnJhcnkodmVnYW4pCmxpYnJhcnkoc2YpCmxpYnJhcnkoY29uY2F2ZW1hbikgI3BvbHlnb24gYXJvdW5kIHBvaW50cwpsaWJyYXJ5KGJldGFwYXJ0KSAjYWxsb3dzIHVzIHRvIHBhcnRpdGlvbiBiZXRhIGRpdmVyc2l0eQpsaWJyYXJ5KGdlb3NwaGVyZSkKbGlicmFyeShnZ3B1YnIpICNzdGF0X3JlZ2xpbmVfZXF1YXRpb24KbGlicmFyeShubG1lKQpsaWJyYXJ5KG1nY3YpICN0byBtYWtlIGdhbQpsaWJyYXJ5KGNvd3Bsb3QpCmxpYnJhcnkobG1lNCkKCiNQdWxsIERpc3NpbWlsYXJpdHkgTWVhbnMKZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yIDwtIHJlYWRSRFMoaGVyZTo6aGVyZSgib3V0cHV0IiwiZGlzc2ltaWxhcml0aWVzIiwgImRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuci5yZHMiKSkKCiNtYWtlIHN1cnZleSBhbmQgc3VydmV5IHVuaXQgZmFjdG9ycwpkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHN1cnZleTo9ZmFjdG9yKHN1cnZleSldWyxzdXJ2ZXlfdW5pdDo9ZmFjdG9yKHN1cnZleV91bml0KV0KCiNhZGp1c3QgeWVhcnMKZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyX2FkaiA6PSB5ZWFyLW1pbih5ZWFyKSsxXQoKI2FkZCBuZXcgdmFyaWFibGUgZm9yIHllYXIgaW4gc2VxdWVuY2UgcGVyIHJlZ2lvbgpkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLGZpcnN0X3llYXIgOj0gbWluKHllYXIpLC4oc3VydmV5X3VuaXQpXQpkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLGxhc3RfeWVhciA6PSBtYXgoeWVhciksLihzdXJ2ZXlfdW5pdCldCgojZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyX2luX3NlcSA6PSB5ZWFyLWZpcnN0X3llYXIrMV0KCmRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcnNfc2FtcGxlZCA6PSBsYXN0X3llYXItZmlyc3RfeWVhcisxXQoKCmBgYAoKIyMjUGFsZXR0ZSBmb3IgUGxvdHRpbmcKUGFsZXR0ZSBmb3IgcGxvdHRpbmcgYWxsIDM3IHN1cnZleSB1bml0cwpgYGB7ciBsaW5rIGNvbG9ycyB0byBzdXJ2ZXkgdW5pdHN9CnN1cnZleV91bml0Lmxpc3QgPC0gbGV2ZWxzKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclssc3VydmV5X3VuaXRdKQoKcGFsZXR0ZV80MiA8LSBjKAogICIjNUE1MTU2IiwgI0FJCiAgIiNERjAwREIiLCAjQklUUy0xCiAgIiNEQjhFREEiLCAjQklUUy00CiAgIiNGNjIyMkUiLCAjQ0hMCiAgIiNGOEExOUYiLCAjREZPLU5GCiAgIiMxNkZGMzIiLCAjREZPLVFDUwogICIjMzI1QTlCIiwgI0VCUwogICIjMzI4M0ZFIiwgI0VWSE9FCiAgIiNGRUFGMTYiLCAjRlItQ0dGUwogICIjZmNjYjZkIiwgI0dNRVgtRmFsbAogICIjMUM4MzU2IiwgI0dNRVgtU3VtbWVyCiAgIiNDNDQ1MUMiLCAjR09BCiAgIiM4NTY2MEQiLCAjR1JMLURFCiAgIiNCMDAwOUYiLCAjR1NMLU4KICAiI0JGNzlCOCIsICNHU0wtUwogICIjMUNCRTRGIiwgI0lDRS1HRlMKICAiIzc4MkFCNiIsICNJRS1JR0ZTCiAgIiM5MEFEMUMiLCAjTUVESVRTCiAgIiM2QjAwM0EiLCAjTkFNCiAgIiNBNzVCMDAiLCAjTkVVUy1GYWxsCiAgIiNFM0IwNzIiLCAjTkVVUy1TcHJpbmcKICAiIzAyRThCNiIsICNOSUdGUy0xCiAgIiM5N0U3RDUiLCAjTklHRlMtNAogICIjQjAwMDY4IiwgI05vci1CVFMtMwogICIjMDBCOUUzIiwgI05TLUlCVFMtMQogICIjOTVFMkY0IiwgI05TLUlCVFMtMwogICIjQjNDRTczIiwgI05aLUNIQVQKICAiIzY4OTUwMCIsICNOWi1FQ1NJCiAgIiMzNjRkMDIiLCNOWi1TVUJBCiAgIiNBQUY0MDAiLCAjTlotV0NTSQogICIjQUEwREZFIiwgI1BULUlCVFMKICAiIzdmOWViOCIsICNST0NLQUxMCiAgIiNGQTAwODciLCAjUy1HRU9SRwogICIjREVBMEZEIiwgI1NDUy1TdW1tZXIKICAiI0ZDRUY4OCIsICNTRVVTLWZhbGwKICAiI0E1OTQwNSIsICNTRVVTLXNwcmluZwogICIjRkNFMTAwIiwgI1NFVVMtc3VtbWVyCiAgIiM1NDQ1NjMiLCAjU1dDLUlCVFMtMQogICIjYTM3ZmM3IiwgI1NXQy1JQlRTLTQKICAiI0MwNzVBNiIsICNXQ0FOTgogICIjQkRDREZGIiwgI1pBRi1BVEwKICAiIzAwM0VGRiIgICNaQUYtSU5ECikKCmNvbG9yX2xpbmsgPC0gZGF0YS50YWJsZShzdXJ2ZXlfdW5pdCA9IHN1cnZleV91bml0Lmxpc3QsaGV4ID0gcGFsZXR0ZV80MikKYGBgCgpBZGQgbmFtZXMgZm9yIHBsb3R0aW5nCmBgYHtyIGFkZCBuYW1lcyBmb3IgcGxvdHRpbmd9CgpuYW1lX2hlbHBlciA8LSBkYXRhLnRhYmxlKFN1cnZleV9OYW1lX1NlYXNvbiA9IGMoIkFsZXV0aWFuIElzbGFuZHMiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiQmFsdGljIFNlYSBRMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJCYWx0aWMgU2VhIFE0IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkNoaWxlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5ld2ZvdW5kbGFuZCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJRdWVlbiBDaGFybG90dGUgU291bmQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRWFzdGVybiBCZXJpbmcgU2VhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkJheSBvZiBCaXNjYXkiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRW5nbGlzaCBDaGFubmVsIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkd1bGYgb2YgTWV4aWNvIFN1bW1lciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHdWxmIG9mIEFsYXNrYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHcmVlbmxhbmQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTiBHdWxmIG9mIFN0LiBMYXdyZW5jZSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTIEd1bGYgb2YgU3QuIExhd3JlbmNlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkljZWxhbmQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSXJpc2ggU2VhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk1lZGl0ZXJyYW5lYW4iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTmFtaWJpYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJORSBVUyBGYWxsIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5FIFVTIFNwcmluZyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOIElyZWxhbmQgUTEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTiBJcmVsYW5kIFE0IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkJhcmVudHMgU2VhIE5vcndheSBRMyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOIFNlYSBRMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOIFNlYSBRMyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJDaGF0aGFtIFJpc2UgTloiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRSBDb2FzdCBTIElzbGFuZCBOWiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJXIENvYXN0IFMgSXNsYW5kIE5aIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlBvcnR1Z2FsIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlMgR2VvcmdpYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU2NvdGlhbiBTaGVsZiBTdW1tZXIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNFIFVTIEZhbGwiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNFIFVTIFNwcmluZyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU0UgVVMgU3VtbWVyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJXIENvYXN0IFVTIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJBdGxhbnRpYyBPY2VhbiBaQSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSW5kaWFuIE9jZWFuIFpBIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiUm9ja2FsbCBQbGF0ZWF1IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTY290bGFuZCBTaGVsZiBTZWEgUTEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNjb3RsYW5kIFNoZWxmIFNlYSBRNCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRmFsa2xhbmQgSXNsYW5kcyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR3VsZiBvZiBNZXhpY28gRmFsbCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiQmFyZW50cyBTZWEgTm9yd2F5IFExIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTdWItQXJjdGljIE5aIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTY290aWFuIFNoZWxmIFNwcmluZyIpLAogICAgICAgICAgICAgICAgICAgICAgICAgIHN1cnZleV91bml0ID0gYygKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJBSSIsICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJCSVRTLTEiLCAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJCSVRTLTQiLCAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJDSEwiLCAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJERk8tTkYiLCAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJERk8tUUNTIiwgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJFQlMiLCAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJFVkhPRSIsICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJGUi1DR0ZTIiwgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHTUVYLVN1bW1lciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR09BIiwgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR1JMLURFIiwgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR1NMLU4iLCAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR1NMLVMiLCAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSUNFLUdGUyIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSUUtSUdGUyIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTUVESVRTIiwgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTkFNIiwgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTkVVUy1GYWxsIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTkVVUy1TcHJpbmciLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5JR0ZTLTEiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5JR0ZTLTQiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5vci1CVFMtMyIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5TLUlCVFMtMSIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5TLUlCVFMtMyIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5aLUNIQVQiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5aLUVDU0kiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5aLVdDU0kiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlBULUlCVFMiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlMtR0VPUkciLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNDUy1TVU1NRVIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNFVVMtZmFsbCIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNFVVMtc3ByaW5nIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTRVVTLXN1bW1lciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiV0NBTk4iLCAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiWkFGLUFUTCIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiWkFGLUlORCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiUk9DS0FMTCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU1dDLUlCVFMtMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU1dDLUlCVFMtNCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRkFMSyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR01FWC1GYWxsIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOb3ItQlRTLTEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5aLVNVQkEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNDUy1TUFJJTkciCiAgICAgICAgICAgICAgICAgICAgICAgICAgKSkKCgpjb2xvcl9saW5rIDwtIGNvbG9yX2xpbmtbbmFtZV9oZWxwZXIsIG9uID0gInN1cnZleV91bml0Il0KCmBgYAoKCgojI01ha2UgR0FNcwoKQnJheSBDdXJ0aXMKYGBge3IgYnJheSBjdXJ0aXMgZ2Ftc30KYnJheV9jdXJ0aXNfYmFsYW5jZWRfZ2FtIDwtIGdhbShicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2JhbGFuY2VkX21lYW4gfiB5ZWFyICsgcyhzdXJ2ZXlfdW5pdCwgeWVhciwgYnMgPSAiZnMiLCBtID0gMSksI3JhbmRvbSBzbW9vdGgKICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRhdGEgPSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIpCgpgYGAKCgojI01ha2UgTE1FUlMKCkJyYXkKKlRoZXNlIGFsbCBjb252ZXJnZWQqCmBgYHtyIGJyYXl9CiNydW5uaW5nIHdpdGggbG1lIGluc3RlYWQgb2YgbG1lciBnYXZlIHNhbWUgcmVzdWx0cywgYnV0IGFsbG93ZWQgZm9yIGNhbGN1bGF0aW9uIG9mIHAtdmFsdWUKYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lIDwtIGxtZShicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2JhbGFuY2VkX21lYW4gfiB5ZWFyX2FkaiwgcmFuZG9tID0gKH4xICsgeWVhcl9hZGp8c3VydmV5X3VuaXQpLGRhdGEgPSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIpCgojYnV0IGFsc28gcnVuIHdpdGggbG1lciBmb3IgY29uZmludApicmF5X2N1cnRpc19iYWxhbmNlZF9sbWVyIDwtIGxtZXIoYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV9iYWxhbmNlZF9tZWFuIH4geWVhcl9hZGogKyAoMSArIHllYXJfYWRqfHN1cnZleV91bml0KSxkYXRhID0gZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yKQoKc3VtbWFyeShicmF5X2N1cnRpc19iYWxhbmNlZF9sbWUpCmFub3ZhKGJyYXlfY3VydGlzX2JhbGFuY2VkX2xtZSkKCmJyYXlfY3VydGlzX2JhbGFuY2VkX2NvZWZzIDwtIGRhdGEudGFibGUoY29lZihicmF5X2N1cnRpc19iYWxhbmNlZF9sbWUpKQpicmF5X2N1cnRpc19iYWxhbmNlZF9jb2Vmc1ssc3VydmV5X3VuaXQgOj0gcm93bmFtZXMoY29lZihicmF5X2N1cnRpc19iYWxhbmNlZF9sbWUpKV1bLFllYXIgOj0gcm91bmQoeWVhcl9hZGosNSldWyxJbnRlcmNlcHQgOj0gcm91bmQoYChJbnRlcmNlcHQpYCwyKV0KVmlldyhicmF5X2N1cnRpc19iYWxhbmNlZF9jb2VmcykKCmJyYXlfY3VydGlzX2JhbGFuY2VkX2NvZWZzIDwtIGJyYXlfY3VydGlzX2JhbGFuY2VkX2NvZWZzW2NvbG9yX2xpbmssIG9uID0gInN1cnZleV91bml0Il0KCmJyYXlfY3VydGlzX2JhbGFuY2VkX2NvZWZzLmV4cCA8LSBicmF5X2N1cnRpc19iYWxhbmNlZF9jb2Vmc1ssLihTdXJ2ZXlfTmFtZV9TZWFzb24sIEludGVyY2VwdCwgWWVhcildCgojZXhwb3J0IHRoaXMgdGFibGUKZndyaXRlKGJyYXlfY3VydGlzX2JhbGFuY2VkX2NvZWZzLmV4cCwgZmlsZSA9IGhlcmU6OmhlcmUoIm91dHB1dCIsImJyYXlfY3VydGlzX2JhbGFuY2VkX2NvZWZzLmV4cC5jc3YiKSkKYGBgCgpHZXQgTE1FUiBtb2RlbCBhcyBwcmVkaWN0aW9ucwpgYGB7cn0KCiMgbmVlZCB0byBzb3J0IG91dCB5ZWFyIGluIHNlcSB2ZXJzdXMgb3ZlcmFsbCB5ZWFyIG1vZGVscwojbmV3IGRhdGEgZm9yIGxtZXIKbG1lcl95ZWFyIDwtIHNlcShtaW4oZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSksIG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSwgYnkgPSAwLjEpCgpsbWVyX3llYXJfYWRqIDwtIHNlcShtaW4oZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyX2Fkal0pLCBtYXgoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyX2Fkal0pLCBieSA9IDAuMSkKCiNwcmVkaWN0IGF2ZXJhZ2UgbG1lcgpsbWVyX2JyYXlfYmFsYW5jZWRfcHJlZGljdGlvbnMgPC0gZGF0YS50YWJsZSh5ZWFyID0gbG1lcl95ZWFyLCB5ZWFyX2FkaiA9IGxtZXJfeWVhcl9hZGopCgojY29uZmlkZW5jZSBpbnRlcnZhbHMKYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcl9jb25maW50IDwtIGNvbmZpbnQoYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcikKCiNwb3B1bGF0ZSBkYXRhIHRhYmxlIG9mIGxtZXIgcHJlZGljdGlvbnMKbG1lcl9icmF5X2JhbGFuY2VkX3ByZWRpY3Rpb25zWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzIDo9IGZpeGVmKGJyYXlfY3VydGlzX2JhbGFuY2VkX2xtZXIpW1sxXV0gKyBmaXhlZihicmF5X2N1cnRpc19iYWxhbmNlZF9sbWVyKVtbMl1dICogeWVhcl9hZGpdWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzX2xvd2VyQ0kgOj0gYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcl9jb25maW50WzVdICsgYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcl9jb25maW50WzZdICogeWVhcl9hZGpdWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzX3VwcGVyQ0kgOj0gYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcl9jb25maW50WzExXSArIGJyYXlfY3VydGlzX2JhbGFuY2VkX2xtZXJfY29uZmludFsxMl0gKiB5ZWFyX2Fkal0KYGBgCgoKCiMjI01vdmUgZm9yd2FyZCB3aXRoIEJyYXkgQ3VydGlzIEJhbGFuY2VkIChvdGhlcnMgdG8gc3VwcGxlbWVudCkKCgpDb2VmZmljaWVudHMgZm9yIExNRVIgYnkgc3VydmV5X3VuaXQKYGBge3J9CiN1bmlxdWUgc3VydmV5IHVuaXQgeWVhciBjb21ib3MKc3VydmV5X3VuaXRfc2FtcGxpbmdfeWVhcnMgPC0gdW5pcXVlKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclssLihzdXJ2ZXlfdW5pdCwgeWVhcl9hZGosIHllYXIsIHllYXJzX3NhbXBsZWQpXSkKCiMgc2VlIGdyb3VwIGNvZWZmaWNpZW50cwptb2RlbF9jb2Vmc19yZWR1Y2VkIDwtIGRhdGEudGFibGUodHJhbnNmb3JtKGFzLmRhdGEuZnJhbWUocmFuZWYoYnJheV9jdXJ0aXNfYmFsYW5jZWRfbG1lcikpLCBsd3IgPSBjb25kdmFsIC0gMS45Nipjb25kc2QsIHVwciA9IGNvbmR2YWwgKyAxLjk2KmNvbmRzZCkpCiNodHRwczovL3N0YWNrb3ZlcmZsb3cuY29tL3F1ZXN0aW9ucy82OTgwNTUzMi9leHRyYWN0LXRoZS1jb25maWRlbmNlLWludGVydmFscy1vZi1sbWVyLXJhbmRvbS1lZmZlY3RzLXBsb3R0ZWQtd2l0aC1kb3RwbG90cmEKCgojT05MWSBTTE9QRVMKbW9kZWxfY29lZnNfcmVkdWNlZCA8LSBtb2RlbF9jb2Vmc19yZWR1Y2VkW3Rlcm0gPT0gInllYXJfYWRqIixdCgptb2RlbF9jb2Vmc19yZWR1Y2VkWyxzdXJ2ZXlfdW5pdCA6PSBncnBdWyx5ZWFyX2FkaiA6PSBjb25kdmFsXQoKI21lcmdlIHdpdGggZHVyYXRpb24gb2Ygc3VydmV5Cm1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoIDwtIG1vZGVsX2NvZWZzX3JlZHVjZWRbc3VydmV5X3VuaXRfc2FtcGxpbmdfeWVhcnMsIG9uID0gInN1cnZleV91bml0Il0KCgptb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aFsseWVhcnNfc2FtcGxlZCA6PSBhcy5udW1lcmljKHllYXJzX3NhbXBsZWQpXVssRGlyZWN0aW9uYWxfQ2hhbmdlIDo9IGlmZWxzZSh5ZWFyX2FkaiA+IDAsICJEaWZmZXJlbnRpYXRpb24iLCJIb21vZ2VuaXphdGlvbiIpXQoKI2RvZXMgaXQgY3Jvc3MgemVybz8KbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGhbLHNpZ25pZmljYW50IDo9IGlmZWxzZShsd3IgPjAgJiB1cHI+MCxULGlmZWxzZShsd3I8MCAmIHVwcjwwLFQsRikpXQoKI3NpZ25pZmljYW50IGRpcmVjdGlvbmFsIGNoYW5nZQptb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aFssRGlyZWN0aW9uYWxfQ2hhbmdlX3NpZyA6PSBpZmVsc2UoeWVhcl9hZGogPiAwICYgc2lnbmlmaWNhbnQgPT0gVCwgIkRpZmZlcmVudGlhdGlvbiIsaWZlbHNlKHllYXJfYWRqIDwgMCAmIHNpZ25pZmljYW50ID09IFQsICJIb21vZ2VuaXphdGlvbiIsICJObyB0cmVuZCBpblxuZGlzc2ltaWxhcml0eSIpKV0KCgojbWluIG1heCBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzCm1pbl9icmF5X3JlZHVjZWQgPC0gbWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclssYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV9iYWxhbmNlZF9tZWFuXSwgbmEucm0gPSBUKQptYXhfYnJheV9yZWR1Y2VkIDwtIG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfYmFsYW5jZWRfbWVhbl0sIG5hLnJtID0gVCkKCm1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoIDwtIG1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoW2NvbG9yX2xpbmssIG9uID0gInN1cnZleV91bml0Il0KCiNkZWxldGUgYW55IE5Bcwptb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCA8LSBuYS5vbWl0KG1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLCBjb2xzID0gInNpZ25pZmljYW50IikKCiNvcmRlciB0YWJsZSBieSBjb2VmZmljaWVudApzZXRvcmRlcihtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCwgeWVhcl9hZGopCgpCQ19iYWxhbmNlZF9tb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aC51bmlxdWUgPC0gdW5pcXVlKG1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoWywuKGNvbmR2YWwsY29uZHNkLCBsd3IsIHVwciwgc3VydmV5X3VuaXQsIHllYXJfYWRqLCB5ZWFyc19zYW1wbGVkLCBEaXJlY3Rpb25hbF9DaGFuZ2UsIGhleCwgU3VydmV5X05hbWVfU2Vhc29uLCBzaWduaWZpY2FudCwgRGlyZWN0aW9uYWxfQ2hhbmdlX3NpZyldKSAKCiNleHRyYWN0IGNvbG9yIGhleGVzCiN5ZWFyIGFkaiBjb2VmIG9yZGVyCmNvbG9yX3llYXJfYWRqX29yZGVyIDwtIEJDX2JhbGFuY2VkX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZVssaGV4XQoKI2FscGhhYmV0aWNhbCBvcmRlcgpCQ19iYWxhbmNlZF9tb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aC51bmlxdWUuYWxwaGEgPC0gc2V0b3JkZXIoQkNfYmFsYW5jZWRfbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgudW5pcXVlLCBTdXJ2ZXlfTmFtZV9TZWFzb24pCgpCQ19iYWxhbmNlZF9tb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aC51bmlxdWUuYWxwaGFbLHRyZW5kX2NvbG9yIDo9IGlmZWxzZShEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnID09ICJIb21vZ2VuaXphdGlvbiIsICIjZTdhYzViIiwgaWZlbHNlKERpcmVjdGlvbmFsX0NoYW5nZV9zaWcgPT0gIkRpZmZlcmVudGlhdGlvbiIsIiM5MWM4NzQiLCIjY2JiZmRlIikpXQoKY29sb3JfYWxwaGFfb3JkZXIgPC0gQkNfYmFsYW5jZWRfbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgudW5pcXVlLmFscGhhWyxoZXhdCmNvbG9yX2FscGhhX29yZGVyX2J5dHJlbmQgPC0gQkNfYmFsYW5jZWRfbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgudW5pcXVlLmFscGhhWywgdHJlbmRfY29sb3JdCgpzYXZlUkRTKEJDX2JhbGFuY2VkX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZSwgaGVyZTo6aGVyZSgib3V0cHV0IiwicmVnaW9uX3N0YXRzIiwiQkNfYmFsYW5jZWRfbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgudW5pcXVlLlJkcyIpKQpgYGAKCkJhciBQbG90IENvZWZmaWNpZW50IExNRVIKYGBge3IgYmFyIHBsb3Qgb2YgY29lZmZpY2llbnRzfQpCQ19CQUxBTkNFRF9EaXNzaW1pbGFyaXR5X0NvZWZfZXJyb3JiYXJfcmVkdWNlZCA8LSBnZ3Bsb3QoKSArCiAgICBnZW9tX2Vycm9yYmFyKGRhdGEgPSBtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCwgYWVzKHggPSByZW9yZGVyKFN1cnZleV9OYW1lX1NlYXNvbiwgeWVhcl9hZGopICwgeSA9IHllYXJfYWRqLCBsYWJlbCA9IFN1cnZleV9OYW1lX1NlYXNvbiwgeW1pbiA9IGx3ciwgeW1heCA9IHVwciksIGZpbGwgPSAiZ3JleSIsIHdpZHRoID0gMCkgKyAjYWRkIGNvbmZpZGVuY2UgaW50ZXJ2YWxzCiAgZ2VvbV9wb2ludChkYXRhID0gbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgsIGFlcyh4ID0gcmVvcmRlcihTdXJ2ZXlfTmFtZV9TZWFzb24sIHllYXJfYWRqKSAsIHkgPSB5ZWFyX2FkaiwgbGFiZWwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24sIHNpemUgPSB5ZWFyc19zYW1wbGVkLCBmaWxsID0gRGlyZWN0aW9uYWxfQ2hhbmdlX3NpZywgY29sb3IgPSBEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnKSwgc3RhdCA9ICdpZGVudGl0eScsIHNoYXBlID0gMjEpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKCJ3aGl0ZSIsImJsYWNrIiwiZ3JleSIpLCBuYW1lID0gIkRpc3NpbWlsYXJpdHkgdHJlbmQiLCBndWlkZT0ibm9uZSIpICsKICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gYygiYmxhY2siLCJibGFjayIsImdyZXkiKSwgbmFtZSA9ICJEaXNzaW1pbGFyaXR5IHRyZW5kIiwgZ3VpZGU9Im5vbmUiKSArCiAgc2NhbGVfc2l6ZV9iaW5uZWQocmFuZ2UgPSBjKDEsOCksIG5hbWUgPSAiU3VydmV5IHBlcmlvZCBsZW5ndGgiKSArCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCkgKwogIHNjYWxlX3lfY29udGludW91cyhicmVha3MgPSBzZXEoLTAuMDA1LCAwLjAwNzUsIGJ5ID0gMC4wMDI1KSwgbGFiZWxzID0gYygiLTAuMDA1IiwiIiwiMCIsICIiLCAiMC4wMDUiLCAgIiIpKSArCiAgeGxhYigiU3VydmV5IHVuaXQiKSArCiAgeWxhYigizrItZGl2ZXJzaXR5IHRyZW5kIikgKyAjYmFsYW5jZWQgQkMgZGlzc2ltaWxhcml0eSB0cmVuZAogIGNvb3JkX2ZsaXAoKSArCiAgdGhlbWVfY2xhc3NpYygpICsKICB0aGVtZShheGlzLnRleHQueSA9IGVsZW1lbnRfdGV4dChjb2xvdXIgPSBjb2xvcl95ZWFyX2Fkal9vcmRlciwgZmFjZSA9ICJib2xkIiksIGF4aXMudGl0bGUueSA9IGVsZW1lbnRfYmxhbmsoKSwgYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGxlZ2VuZC5wb3NpdGlvbiA9IGMoMC4yNSwwLjgpLCBsZWdlbmQuZGlyZWN0aW9uID0gInZlcnRpY2FsIikKCiNwdWxsIGxlZ2VuZCBmb3IgaG9tb2dlbml6YXRpb24KZGlyZWN0aW9uYWxfY2hhbmdlX2xlZ2VuZF9wbG90IDwtIEJDX0JBTEFOQ0VEX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkICsgCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYygid2hpdGUiLCJibGFjayIsImdyZXkiKSwgbmFtZSA9ICJEaXNzaW1pbGFyaXR5IHRyZW5kIikgKwogIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSBjKCJibGFjayIsImJsYWNrIiwiZ3JleSIpLCBuYW1lID0gIkRpc3NpbWlsYXJpdHkgdHJlbmQiKSArCiAgc2NhbGVfc2l6ZV9iaW5uZWQocmFuZ2UgPSBjKDEsOCksIG5hbWUgPSAiU3VydmV5IHBlcmlvZCBsZW5ndGgiLCBndWlkZSA9ICJub25lIikgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJyaWdodCIsIGxlZ2VuZC5iYWNrZ3JvdW5kID0gZWxlbWVudF9yZWN0KGZpbGw9ICJ0cmFuc3BhcmVudCIpKSArCiAgZ3VpZGVzKGNvbG91ciA9IGd1aWRlX2xlZ2VuZChvdmVycmlkZS5hZXMgPSBsaXN0KHNpemU9NikpKQoKCkJDX0JBTEFOQ0VEX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkCgojQUxUIGdyZXkgc2NhbGUKQkNfQkFMQU5DRURfRGlzc2ltaWxhcml0eV9Db2VmX2Vycm9yYmFyX3JlZHVjZWRfZ3JleXNjYWxlIDwtIGdncGxvdCgpICsKICAgIGdlb21fZXJyb3JiYXIoZGF0YSA9IG1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLCBhZXMoeCA9IHJlb3JkZXIoU3VydmV5X05hbWVfU2Vhc29uLCB5ZWFyX2FkaikgLCB5ID0geWVhcl9hZGosIGxhYmVsID0gU3VydmV5X05hbWVfU2Vhc29uLCB5bWluID0gbHdyLCB5bWF4ID0gdXByKSwgZmlsbCA9ICJncmV5Iiwgd2lkdGggPSAwKSArICNhZGQgY29uZmlkZW5jZSBpbnRlcnZhbHMKICBnZW9tX3BvaW50KGRhdGEgPSBtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCwgYWVzKHggPSByZW9yZGVyKFN1cnZleV9OYW1lX1NlYXNvbiwgeWVhcl9hZGopICwgeSA9IHllYXJfYWRqLCBsYWJlbCA9IFN1cnZleV9OYW1lX1NlYXNvbiwgc2l6ZSA9IHllYXJzX3NhbXBsZWQsIGZpbGwgPSBEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnLCBjb2xvciA9IERpcmVjdGlvbmFsX0NoYW5nZV9zaWcpLCBzdGF0ID0gJ2lkZW50aXR5Jywgc2hhcGUgPSAyMSkgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGMoIndoaXRlIiwiYmxhY2siLCJncmV5IiksIG5hbWUgPSAiRGlzc2ltaWxhcml0eSB0cmVuZCIsIGd1aWRlPSJub25lIikgKwogIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSBjKCJibGFjayIsImJsYWNrIiwiZ3JleSIpLCBuYW1lID0gIkRpc3NpbWlsYXJpdHkgdHJlbmQiLCBndWlkZT0ibm9uZSIpICsKICBzY2FsZV9zaXplX2Jpbm5lZChyYW5nZSA9IGMoMSw4KSwgbmFtZSA9ICJTdXJ2ZXkgcGVyaW9kIGxlbmd0aCIpICsKICBnZW9tX2hsaW5lKHlpbnRlcmNlcHQgPSAwKSArCiAgc2NhbGVfeV9jb250aW51b3VzKGJyZWFrcyA9IHNlcSgtMC4wMDUsIDAuMDA3NSwgYnkgPSAwLjAwMjUpLCBsYWJlbHMgPSBjKCItMC4wMDUiLCIiLCIwIiwgIiIsICIwLjAwNSIsICAiIikpICsKICB4bGFiKCJTdXJ2ZXkgdW5pdCIpICsKICB5bGFiKCLOsi1kaXZlcnNpdHkgdHJlbmQiKSArICNiYWxhbmNlZCBCQyBkaXNzaW1pbGFyaXR5IHRyZW5kCiAgY29vcmRfZmxpcCgpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIHRoZW1lKGF4aXMudGV4dC55ID0gZWxlbWVudF90ZXh0KGZhY2UgPSAiYm9sZCIpLCBheGlzLnRpdGxlLnkgPSBlbGVtZW50X2JsYW5rKCksIGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBsZWdlbmQucG9zaXRpb24gPSBjKDAuMjUsMC44KSwgbGVnZW5kLmRpcmVjdGlvbiA9ICJ2ZXJ0aWNhbCIpCmBgYApBbHRlcm5hdGl2ZWx5LCB3ZSBjb2xvciAgdGhpcyBwbG90IGJ5IHRyZW5kIGV4cGVyaWVuY2VkCgpgYGB7cn0KIyIjNzNCQTREIiwiI0UwOTYyQyIsIiNjYmJmZGUiCgpCQ19iYWxhbmNlZF9EaXNzaW1pbGFyaXR5X0NvZWZfZXJyb3JiYXJfcmVkdWNlZF9jb2xvcmJ5dHJlbmQgPC0gZ2dwbG90KCkgKwogICAgZ2VvbV9lcnJvcmJhcihkYXRhID0gbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgsIGFlcyh4ID0gcmVvcmRlcihTdXJ2ZXlfTmFtZV9TZWFzb24sIHllYXJfYWRqKSAsIHkgPSB5ZWFyX2FkaiwgbGFiZWwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24sIHltaW4gPSBsd3IsIHltYXggPSB1cHIpLCBmaWxsID0gImdyZXkiLCB3aWR0aCA9IDApICsgI2FkZCBjb25maWRlbmNlIGludGVydmFscwogIGdlb21fcG9pbnQoZGF0YSA9IG1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLCBhZXMoeCA9IHJlb3JkZXIoU3VydmV5X05hbWVfU2Vhc29uLCB5ZWFyX2FkaikgLCB5ID0geWVhcl9hZGosIGxhYmVsID0gU3VydmV5X05hbWVfU2Vhc29uLCBzaXplID0geWVhcnNfc2FtcGxlZCwgY29sb3IgPSBEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnKSwgc3RhdCA9ICdpZGVudGl0eScpICsKICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gYygiIzczQkE0RCIsIiNFMDk2MkMiLCIjY2JiZmRlIiksIG5hbWUgPSAiRGlzc2ltaWxhcml0eSB0cmVuZCIsIGd1aWRlPSJub25lIikgKwogIHNjYWxlX3NpemVfYmlubmVkKHJhbmdlID0gYygxLDgpLCBuYW1lID0gIlN1cnZleSBwZXJpb2QgbGVuZ3RoXG4oeWVhcnMpIikgKwogIGdlb21faGxpbmUoeWludGVyY2VwdCA9IDApICsKICBzY2FsZV95X2NvbnRpbnVvdXMoYnJlYWtzID0gc2VxKC0wLjAwNSwgMC4wMDc1LCBieSA9IDAuMDAyNSksIGxhYmVscyA9IGMoIi0wLjAwNSIsIiIsIjAiLCAiIiwgIjAuMDA1IiwgICIiKSkgKwogIHhsYWIoIlN1cnZleSB1bml0IikgKwogIHlsYWIoIs6yLWRpdmVyc2l0eSB0cmVuZCIpICsgI2JhbGFuY2VkIEJDIGRpc3NpbWlsYXJpdHkgdHJlbmQKICBjb29yZF9mbGlwKCkgKwogIHRoZW1lX2NsYXNzaWMoKSArCiAgdGhlbWUoYXhpcy50ZXh0LnkgPSBlbGVtZW50X3RleHQoZmFjZSA9ICJib2xkIiksIGF4aXMudGl0bGUueSA9IGVsZW1lbnRfYmxhbmsoKSwgYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGxlZ2VuZC5wb3NpdGlvbiA9IGMoMC4zLDAuOCksIGxlZ2VuZC5kaXJlY3Rpb24gPSAidmVydGljYWwiLCBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE2KSkKCmRpcmVjdGlvbmFsX2NoYW5nZV9sZWdlbmRfcGxvdF9jb2xvcmJ5dHJlbmQgPC0gQkNfYmFsYW5jZWRfRGlzc2ltaWxhcml0eV9Db2VmX2Vycm9yYmFyX3JlZHVjZWRfY29sb3JieXRyZW5kICsgCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gInJpZ2h0IiwgbGVnZW5kLmJhY2tncm91bmQgPSBlbGVtZW50X3JlY3QoZmlsbD0gInRyYW5zcGFyZW50IiksIAogICAgICAgICBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE2KSkgKwogIGd1aWRlcyhjb2xvdXIgPSBndWlkZV9sZWdlbmQob3ZlcnJpZGUuYWVzID0gbGlzdChzaXplPTYpKSwgc2l6ZSA9ICJub25lIikKYGBgCgpXYXZ5IExpbmUgUGxvdCBmb3IgR0FNcwoKR2VuZXJhdGUgcHJlZGljdGVkIHZhbHVlcwpgYGB7ciBnZW5lcmF0ZSBwcmVkaWN0ZWQgdmFsdWVzIEdBTX0KCiNhZGQgY29sb3JzIGFuZCBuYW1lcyB0byBmdWxsIGRpc3NpbWlsYXJpdHkgZGF0YSB0YWJsZQpkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIgPC0gZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yW2NvbG9yX2xpbmssIG9uID0gInN1cnZleV91bml0Il0KCiNnZW5lcmF0ZSBuZXcgZGF0YSB0byBzbW9vdGggbGluZXMgKG5lZWQgeWVhciBhbmQgc2Vhc29uIHN1cnZleSBjb21iaW5hdGlvbnMpCnllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0IDwtIGRhdGEudGFibGUoc3VydmV5X3VuaXQgPSBhcy5jaGFyYWN0ZXIoTlVMTCksIHllYXIgPSBhcy5udW1lcmljKE5VTEwpLCB5ZWFyX2FkaiA9IGFzLm51bWVyaWMoTlVMTCApKQoKZm9yIChpIGluIDE6bGVuZ3RoKHN1cnZleV91bml0Lmxpc3QpKSB7CiAgI2dlbmVyYXRlIHllYXIgdmVjdG9ycwogIHN1cnZleV91bml0X3llYXJzIDwtIHVuaXF1ZShkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbc3VydmV5X3VuaXQgPT0gc3VydmV5X3VuaXQubGlzdFtpXSwuKHN1cnZleV91bml0LCB5ZWFyLCB5ZWFyX2FkaildKQogIAogIHllYXJzIDwtIHNlcShtaW4oc3VydmV5X3VuaXRfeWVhcnNbLHllYXJdKSwgbWF4KHN1cnZleV91bml0X3llYXJzWyx5ZWFyXSksIGJ5ID0gMC4xKQogIAogIHllYXJfYWRqdXN0IDwtIHNlcShtaW4oc3VydmV5X3VuaXRfeWVhcnNbLHllYXJfYWRqXSksIG1heChzdXJ2ZXlfdW5pdF95ZWFyc1sseWVhcl9hZGpdKSwgYnkgPSAwLjEpCiAgCiAgeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRfYWRkaXRpb24gPC0gZGF0YS50YWJsZShzdXJ2ZXlfdW5pdCA9IHN1cnZleV91bml0Lmxpc3RbaV0sIHllYXIgPSB5ZWFycywgeWVhcl9hZGogPSB5ZWFyX2FkanVzdCkKICAKICB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCA8LSByYmluZCh5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCwgeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRfYWRkaXRpb24pCn0KCiNhZGQgY29sb3JzIGFuZCBuYW1lcyB0byBmdWxsIHllYXIgYW5kIHN1cnZleSB1bml0IGNvbWJpbmF0aW9uIHRhYmxlCnllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0IDwtIHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0W2NvbG9yX2xpbmssIG9uID0gInN1cnZleV91bml0Il0KYGBgCgoKCkdldCBtb2RlbCBhcyBwcmVkaWN0aW9ucwpgYGB7cn0KI2ZvciBwbG90dGluZywgZ2V0IG1vZGVsIGFzIHByZWRpY3Rpb25zCmJyYXlfY3VydGlzX2JhbGFuY2VkX2dhbV9wcmVkaWN0aW9ucyA8LSBwcmVkaWN0KGJyYXlfY3VydGlzX2JhbGFuY2VkX2dhbSwgc2UuZml0ID0gVFJVRSwgbmV3ZGF0YSA9IHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0KQoKI21lcmdlIGludG8gdGFibGUKeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRbLGJyYXlfZ2xtX21vZF9maXQgOj0gYnJheV9jdXJ0aXNfYmFsYW5jZWRfZ2FtX3ByZWRpY3Rpb25zJGZpdF1bLGJyYXlfZ2xtX21vZF9maXRfU0UgOj0gYnJheV9jdXJ0aXNfYmFsYW5jZWRfZ2FtX3ByZWRpY3Rpb25zJHNlLmZpdF0KCmBgYAoKClByb2R1Y2UgUGxvdCBvZiBHQU0gYW5kIG1lYW4gTE1FUiBsaW5lCmBgYHtyIHBsb3QgR0FNIGFuZCBtZWFuIExNRVIgbGluZXN9CnBvaW50c193YXZ5bGluZXNfYnJheV9iYWxhbmNlZF95ZWFyX3JlZHVjZWRfZ2FtX25vbG1lciA8LSBnZ3Bsb3QoKSArCiAjIGdlb21fcmliYm9uKGRhdGEgPSBsbWVyX2JyYXlfYmFsYW5jZWRfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeW1pbiA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfbG93ZXJDSSwgeW1heCA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfdXBwZXJDSSksIGZpbGwgPSAiZ3JleSIsIGFscGhhID0gMC4yKSArCiAgZ2VvbV9wb2ludChkYXRhID0gbmEub21pdChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIsY29scyA9ICJ5ZWFyX2FkaiIpLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2JhbGFuY2VkX21lYW4sCiAgICAgICAgICAgICAgICAgY29sb3IgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYSA9IDAuNSwgc2l6ZSA9IDEpICsKICAgIGdlb21fbGluZShkYXRhID0gbmEub21pdCh5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCxjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQsCiAgICAgICAgICAgICAgICAgY29sb3IgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pKSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IG5hLm9taXQoeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQsY29scyA9ICJ5ZWFyX2FkaiIpLCBhZXMoeCA9IHllYXIsIHltaW49YnJheV9nbG1fbW9kX2ZpdC1icmF5X2dsbV9tb2RfZml0X1NFLCB5bWF4PWJyYXlfZ2xtX21vZF9maXQrYnJheV9nbG1fbW9kX2ZpdF9TRSwgZmlsbCA9ICBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYT0wLjEpICsgI2FkZCBzdGFuZGFyZCBlcnJvcgogIyBnZW9tX2xpbmUoZGF0YSA9IGxtZXJfYnJheV9iYWxhbmNlZF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkcyksIGNvbG9yID0gImJsYWNrIikgKwogICAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9ICBjb2xvcl9hbHBoYV9vcmRlciwgbmFtZSA9ICJTdXJ2ZXkgVW5pdCIpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSAgY29sb3JfYWxwaGFfb3JkZXIsIGd1aWRlID0gIm5vbmUiKSArCiAgdGhlbWVfY2xhc3NpYygpICsKICBsaW1zKHggPSBjKG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSxtYXgoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSkpLAogICAgICAgeSA9IGMoMC4xNSwwLjkpKSArCiAgeGxhYigiWWVhciIpICsKeWxhYigiQmFsYW5jZWQgQkMgZGlzc2ltaWxhcml0eSIpICsKICB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibnVsbCIpCgpwb2ludHNfd2F2eWxpbmVzX2JyYXlfYmFsYW5jZWRfeWVhcl9yZWR1Y2VkX2dhbV9ub2xtZXIKCmdnc2F2ZShwb2ludHNfd2F2eWxpbmVzX2JyYXlfYmFsYW5jZWRfeWVhcl9yZWR1Y2VkX2dhbV9ub2xtZXIsIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lID0icG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW1fbm9sbWVyLmpwZyIsIGhlaWdodCA9IDUsIHdpZHRoID0gNSwgdW5pdCA9ICJpbiIpCgojd2l0aCBsbWVyCgpwb2ludHNfd2F2eWxpbmVzX2JyYXlfYmFsYW5jZWRfeWVhcl9yZWR1Y2VkX2dhbSA8LSBnZ3Bsb3QoKSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IGxtZXJfYnJheV9iYWxhbmNlZF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5bWluID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc19sb3dlckNJLCB5bWF4ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc191cHBlckNJKSwgZmlsbCA9ICJncmV5IiwgYWxwaGEgPSAwLjMpICsKICBnZW9tX3BvaW50KGRhdGEgPSBuYS5vbWl0KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuciwgY29scyA9ICJ5ZWFyX2FkaiIpLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2JhbGFuY2VkX21lYW4sCiAgICAgICAgICAgICAgICAgZmlsbCA9IFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhID0gMC40LCBzaXplID0gMSwgc2hhcGUgPSAyMSwgY29sb3IgPSAid2hpdGUiKSArCiAgICBnZW9tX2xpbmUoZGF0YSA9IG5hLm9taXQoeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQsIGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9nbG1fbW9kX2ZpdCwKICAgICAgICAgICAgICAgICBjb2xvciA9IFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhID0gMC42KSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IG5hLm9taXQoeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQsIGNvbHMgPSAieWVhcl9hZGoiKSwgYWVzKHggPSB5ZWFyLCB5bWluPWJyYXlfZ2xtX21vZF9maXQtYnJheV9nbG1fbW9kX2ZpdF9TRSwgeW1heD1icmF5X2dsbV9tb2RfZml0K2JyYXlfZ2xtX21vZF9maXRfU0UsIGZpbGwgPSAgU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGE9MC4xKSArICNhZGQgc3RhbmRhcmQgZXJyb3IKICBnZW9tX2xpbmUoZGF0YSA9IGxtZXJfYnJheV9iYWxhbmNlZF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkcyksIGNvbG9yID0gImJsYWNrIikgKwogICAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9ICBjb2xvcl9hbHBoYV9vcmRlciwgbmFtZSA9ICJTdXJ2ZXkgVW5pdCIpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSAgY29sb3JfYWxwaGFfb3JkZXIsIGd1aWRlID0gIm5vbmUiKSArCiAgdGhlbWVfY2xhc3NpYygpICsKICBsaW1zKHggPSBjKG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSxtYXgoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSkpLAogICAgICAgeSA9IGMoMC4xNSwwLjkpKSArCiAgeGxhYigiWWVhciIpICsKeWxhYigizrItZGl2ZXJzaXR5IikgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJudWxsIiwgYXhpcy50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGF4aXMudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSkKCnBvaW50c193YXZ5bGluZXNfYnJheV9iYWxhbmNlZF95ZWFyX3JlZHVjZWRfZ2FtCgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW0sIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lID0icG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW0uanBnIiwgaGVpZ2h0ID0gNiwgd2lkdGggPSA2LCB1bml0ID0gImluIikKCiNBTFQKI3Bsb3QgYWxsLCBidXQgc2FtZSBjb2xvciBzY2hlbWUgKGdyZXkpCnBvaW50c193YXZ5bGluZXNfYnJheV9iYWxhbmNlZF95ZWFyX3JlZHVjZWRfZ2FtX2dyZXlzY2FsZSA8LSBnZ3Bsb3QoKSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IGxtZXJfYnJheV9iYWxhbmNlZF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5bWluID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc19sb3dlckNJLCB5bWF4ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc191cHBlckNJKSwgZmlsbCA9ICJncmV5IiwgYWxwaGEgPSAwLjMpICsKICBnZW9tX3BvaW50KGRhdGEgPSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIsCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfYmFsYW5jZWRfbWVhbiwKICAgICAgICAgICAgICAgICBmaWxsID0gU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGEgPSAwLjQsIHNpemUgPSAxLCBzaGFwZSA9IDIxLCBjb2xvciA9ICJ3aGl0ZSIpICsKICAgIGdlb21fbGluZShkYXRhID0geWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQsCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQsCiAgICAgICAgICAgICAgICAgY29sb3IgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYSA9IDAuNikgKwogIGdlb21fcmliYm9uKGRhdGEgPSB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCwgYWVzKHggPSB5ZWFyLCB5bWluPWJyYXlfZ2xtX21vZF9maXQtYnJheV9nbG1fbW9kX2ZpdF9TRSwgeW1heD1icmF5X2dsbV9tb2RfZml0K2JyYXlfZ2xtX21vZF9maXRfU0UsIGZpbGwgPSAgU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGE9MC4xKSArICNhZGQgc3RhbmRhcmQgZXJyb3IKICBnZW9tX2xpbmUoZGF0YSA9IGxtZXJfYnJheV9iYWxhbmNlZF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkcyksIGNvbG9yID0gImJsYWNrIikgKwogICAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9ICByZXAoImJsYWNrIiwgdGltZXMgPSBsZW5ndGgodW5pcXVlKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuciRTdXJ2ZXlfTmFtZV9TZWFzb24pKSksIG5hbWUgPSAiU3VydmV5IFVuaXQiKSArCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gIHJlcCgiYmxhY2siLCB0aW1lcyA9IGxlbmd0aCh1bmlxdWUoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yJFN1cnZleV9OYW1lX1NlYXNvbikpKSwgZ3VpZGUgPSAibm9uZSIpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGxpbXMoeCA9IGMobWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSksCiAgICAgICB5ID0gYygwLjE1LDAuOSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCLOsi1kaXZlcnNpdHkiKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiLCBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKQoKcG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW1fZ3JleXNjYWxlCgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW1fZ3JleXNjYWxlLCBwYXRoID0gaGVyZTo6aGVyZSgiZmlndXJlcyIpLCBmaWxlbmFtZSA9InBvaW50c193YXZ5bGluZXNfYnJheV9iYWxhbmNlZF95ZWFyX3JlZHVjZWRfZ2FtX2dyZXlzY2FsZS5qcGciLCBoZWlnaHQgPSA2LCB3aWR0aCA9IDYsIHVuaXQgPSAiaW4iKQoKYGBgCgpBbHRlcm5hdGl2ZSwgY29sb3IgYnkgdHJlbmQKYGBge3IgY29sb3Igd2F2eSBsaW5lcyBieSB0cmVuZH0KCnBvaW50c193YXZ5bGluZXNfYnJheV9iYWxhbmNlZF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZCA8LSBnZ3Bsb3QoKSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IGxtZXJfYnJheV9iYWxhbmNlZF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5bWluID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc19sb3dlckNJLCB5bWF4ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc191cHBlckNJKSwgZmlsbCA9ICJncmV5IiwgYWxwaGEgPSAwLjMpICsKICBnZW9tX3BvaW50KGRhdGEgPSBuYS5vbWl0KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuciwgY29scyA9ICJ5ZWFyX2FkaiIpLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2JhbGFuY2VkX21lYW4sCiAgICAgICAgICAgICAgICAgZmlsbCA9IFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhID0gMC40LCBzaXplID0gMSwgc2hhcGUgPSAyMSwgY29sb3IgPSAid2hpdGUiKSArCiAgICBnZW9tX2xpbmUoZGF0YSA9IG5hLm9taXQoeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQsIGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9nbG1fbW9kX2ZpdCwKICAgICAgICAgICAgICAgICBjb2xvciA9IFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhID0gMC42KSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IG5hLm9taXQoeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQsIGNvbHMgPSAieWVhcl9hZGoiKSwgYWVzKHggPSB5ZWFyLCB5bWluPWJyYXlfZ2xtX21vZF9maXQtYnJheV9nbG1fbW9kX2ZpdF9TRSwgeW1heD1icmF5X2dsbV9tb2RfZml0K2JyYXlfZ2xtX21vZF9maXRfU0UsIGZpbGwgPSAgU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGE9MC4xKSArICNhZGQgc3RhbmRhcmQgZXJyb3IKICBnZW9tX2xpbmUoZGF0YSA9IGxtZXJfYnJheV9iYWxhbmNlZF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkcyksIGNvbG9yID0gImJsYWNrIikgKwogICAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9ICBjb2xvcl9hbHBoYV9vcmRlcl9ieXRyZW5kLCBuYW1lID0gIlN1cnZleSBVbml0IikgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9ICBjb2xvcl9hbHBoYV9vcmRlcl9ieXRyZW5kLCBndWlkZSA9ICJub25lIikgKwogIHRoZW1lX2NsYXNzaWMoKSArCiAgbGltcyh4ID0gYyhtaW4oZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSksbWF4KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pKSkgKwogIHhsYWIoIlllYXIiKSArCnlsYWIoIs6yLWRpdmVyc2l0eSIpICsKICB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibnVsbCIsIGF4aXMudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBheGlzLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSkpCgpwb2ludHNfd2F2eWxpbmVzX2JyYXlfYmFsYW5jZWRfeWVhcl9yZWR1Y2VkX2dhbV9jb2xvcmJ5dHJlbmQKCmdnc2F2ZShwb2ludHNfd2F2eWxpbmVzX2JyYXlfYmFsYW5jZWRfeWVhcl9yZWR1Y2VkX2dhbV9jb2xvcmJ5dHJlbmQsIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lID0icG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kLmpwZyIsIGhlaWdodCA9IDYsIHdpZHRoID0gNiwgdW5pdCA9ICJpbiIpCmBgYAoKYGBge3J9CiNwbG90IGVhY2ggaW5kZXBlbmRlbnRseSBmb3Igc3VwcGxlbWVudAojYWxsIHN1cnZleSBuYW1lcyA9IAphbGxfc3VydmV5X25hbWVzIDwtIHNvcnQodW5pcXVlKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuciRTdXJ2ZXlfTmFtZV9TZWFzb24pKQojbGlzdCBvZiBwbG90cwpwb2ludHNfd2F2eWxpbmVzX2JyYXlfYmFsYW5jZWRfeWVhcl9yZWR1Y2VkX2dhbV9pbmRpdmlkdWFsIDwtIGxpc3QoKQpmb3IgKGkgaW4gMTpsZW5ndGgoYWxsX3N1cnZleV9uYW1lcykpIHsKcG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbFtbaV1dIDwtIGdncGxvdCgpICsKICBnZW9tX3BvaW50KGRhdGEgPSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbU3VydmV5X05hbWVfU2Vhc29uID09IGFsbF9zdXJ2ZXlfbmFtZXNbaV1dLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X2JhbGFuY2VkX21lYW4pLCBhbHBoYSA9IDAuNCwgY29sb3IgPSAiYmxhY2siKSArCiAgICBnZW9tX2xpbmUoZGF0YSA9IHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0W1N1cnZleV9OYW1lX1NlYXNvbiA9PSBhbGxfc3VydmV5X25hbWVzW2ldXSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9nbG1fbW9kX2ZpdCksIGFscGhhID0gMC42KSArCiAgZ2VvbV9yaWJib24oZGF0YSA9IHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0W1N1cnZleV9OYW1lX1NlYXNvbiA9PSBhbGxfc3VydmV5X25hbWVzW2ldXSwgYWVzKHggPSB5ZWFyLCB5bWluPWJyYXlfZ2xtX21vZF9maXQtYnJheV9nbG1fbW9kX2ZpdF9TRSwgeW1heD1icmF5X2dsbV9tb2RfZml0K2JyYXlfZ2xtX21vZF9maXRfU0UpLCBhbHBoYT0wLjEpICsgI2FkZCBzdGFuZGFyZCBlcnJvcgogIHRoZW1lX2NsYXNzaWMoKSArCiMgIGxpbXMoeCA9IGMobWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucltTdXJ2ZXlfTmFtZV9TZWFzb24gPT0gYWxsX3N1cnZleV9uYW1lc1tpXSx5ZWFyXSksbWF4KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucltTdXJ2ZXlfTmFtZV9TZWFzb24gPT0gYWxsX3N1cnZleV9uYW1lc1tpXSx5ZWFyXSkpLAojICAgICAgIHkgPSBjKDAuMTUsMC45KSkgKwogIHhsYWIoIlllYXIiKSArCnlsYWIoImJldGEtZGl2ZXJzaXR5IikgKwogIGZhY2V0X3dyYXAoflN1cnZleV9OYW1lX1NlYXNvbiwgbmNvbCA9IDUpICsKICB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibnVsbCIsIGF4aXMudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBheGlzLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSkpCgpwcmludChwb2ludHNfd2F2eWxpbmVzX2JyYXlfYmFsYW5jZWRfeWVhcl9yZWR1Y2VkX2dhbV9pbmRpdmlkdWFsW1tpXV0pCgp9CnNhdmVSRFMocG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbCwgaGVyZTo6aGVyZSgiZmlndXJlcyIsInBvaW50c193YXZ5bGluZXNfYnJheV9iYWxhbmNlZF95ZWFyX3JlZHVjZWRfZ2FtX2luZGl2aWR1YWwuUmRzIikpCgojcHJpbnQgdG8gcGRmCmxpYnJhcnkoZ3JpZEV4dHJhKQoKZ2dzYXZlKAogICBmaWxlbmFtZSA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiLCJwb2ludHNfd2F2eWxpbmVzX2JyYXlfYmFsYW5jZWRfeWVhcl9yZWR1Y2VkX2dhbV9pbmRpdmlkdWFsLnBkZiIpLCAKICAgcGxvdCA9IG1hcnJhbmdlR3JvYihwb2ludHNfd2F2eWxpbmVzX2JyYXlfYmFsYW5jZWRfeWVhcl9yZWR1Y2VkX2dhbV9pbmRpdmlkdWFsLCBucm93PTEsIG5jb2w9MSksIAogICB3aWR0aCA9IDguNSwgaGVpZ2h0ID0gMTEKKQoKCmBgYAoKTWVyZ2UgQkMgdmVyc3VzIFllYXIgcGxvdCB3aXRoIEdBTVMgYW5kIFJlZ2lvbiB2cy4gY29lZmZpY2llbnQgcGxvdCBmb3IgTE1FUnMKCmBgYHtyfQoKQkNfQkFMQU5DRURfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kIDwtIGdnZHJhdyh4bGltID0gYygwLCA0MC41KSwgeWxpbSA9IGMoMCwgMjEpKSArCiAgICBkcmF3X3Bsb3QocG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW0sCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeCA9IDEsIHkgPSAxLCB3aWR0aCA9IDIwLCBoZWlnaHQgPSAyMCkgKwogICAgZHJhd19wbG90KEJDX0JBTEFOQ0VEX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkICsKICAgICAgICB0aGVtZShsZWdlbmQua2V5LnNpemUgPSB1bml0KDAuNSwgJ2NtJyksICNjaGFuZ2UgbGVnZW5kIGtleSBzaXplCiAgICAgICAgbGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTYpLCAjY2hhbmdlIGxlZ2VuZCB0aXRsZSBmb250IHNpemUKICAgICAgICBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplPTE0KSksICNjaGFuZ2UgbGVnZW5kIHRleHQgZm9udCBzaXplKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gMjAsIHkgPSAxLCB3aWR0aCA9MTksIGhlaWdodCA9IDIwKSArCiAgICBkcmF3X3Bsb3QoZ2V0X2xlZ2VuZChkaXJlY3Rpb25hbF9jaGFuZ2VfbGVnZW5kX3Bsb3QgKyAKICAgICAgdGhlbWUobGVnZW5kLmtleS5zaXplID0gdW5pdCgwLjUsICdjbScpLCAjY2hhbmdlIGxlZ2VuZCBrZXkgc2l6ZQogICAgICAgIGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplPTE1KSwgI2NoYW5nZSBsZWdlbmQgdGl0bGUgZm9udCBzaXplCiAgICAgICAgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZT0xMykpKSwgI2NoYW5nZSBsZWdlbmQgdGV4dCBmb250IHNpemUpCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAyNywgeSA9IDEyLCB3aWR0aCA9IDMsIGhlaWdodCA9IDIpICsKICBnZW9tX3RleHQoYWVzKHggPSAyLCB5ID0gMjAuNyksIGxhYmVsID0gKCJhLiIpLCBzaXplID04LCBmb250ZmFjZSA9ICJib2xkIikgKwogIGdlb21fdGV4dChhZXMoeCA9MjAsIHkgPSAyMC43KSwgbGFiZWwgPSAoImIuIiksIHNpemUgPTgsIGZvbnRmYWNlID0gImJvbGQiKQoKCmdnc2F2ZShCQ19CQUxBTkNFRF9HQU1fTE1FUl9tZXJnZV9sZWdlbmQsIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lID0gIkJDX0JBTEFOQ0VEX0dBTV9MTUVSX21lcmdlX2xlZ2VuZC5wbmciLCBoZWlnaHQgPSA4LCB3aWR0aCA9IDE0LCB1bml0cyA9ICJpbiIpCgojQUxUIEdSRVkgU0NBTEUKQkNfQkFMQU5DRURfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kX2dyZXlzY2FsZSA8LSBnZ2RyYXcoeGxpbSA9IGMoMCwgNDAuNSksIHlsaW0gPSBjKDAsIDIxKSkgKwogICAgZHJhd19wbG90KHBvaW50c193YXZ5bGluZXNfYnJheV9iYWxhbmNlZF95ZWFyX3JlZHVjZWRfZ2FtX2dyZXlzY2FsZSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gMSwgeSA9IDEsIHdpZHRoID0gMjAsIGhlaWdodCA9IDIwKSArCiAgICBkcmF3X3Bsb3QoQkNfQkFMQU5DRURfRGlzc2ltaWxhcml0eV9Db2VmX2Vycm9yYmFyX3JlZHVjZWRfZ3JleXNjYWxlICsKICAgICAgICB0aGVtZShsZWdlbmQua2V5LnNpemUgPSB1bml0KDAuNSwgJ2NtJyksICNjaGFuZ2UgbGVnZW5kIGtleSBzaXplCiAgICAgICAgbGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTYpLCAjY2hhbmdlIGxlZ2VuZCB0aXRsZSBmb250IHNpemUKICAgICAgICBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplPTE0KSksICNjaGFuZ2UgbGVnZW5kIHRleHQgZm9udCBzaXplKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gMjAsIHkgPSAxLCB3aWR0aCA9IDE5LCBoZWlnaHQgPSAyMCkgKwogICAgZHJhd19wbG90KGdldF9sZWdlbmQoZGlyZWN0aW9uYWxfY2hhbmdlX2xlZ2VuZF9wbG90ICsgCiAgICAgIHRoZW1lKGxlZ2VuZC5rZXkuc2l6ZSA9IHVuaXQoMC41LCAnY20nKSwgI2NoYW5nZSBsZWdlbmQga2V5IHNpemUKICAgICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xNSksICNjaGFuZ2UgbGVnZW5kIHRpdGxlIGZvbnQgc2l6ZQogICAgICAgIGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemU9MTMpKSksICNjaGFuZ2UgbGVnZW5kIHRleHQgZm9udCBzaXplKQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAyNywgeSA9IDEyLCB3aWR0aCA9IDMsIGhlaWdodCA9IDIpICsKICBnZW9tX3RleHQoYWVzKHggPSAyLCB5ID0gMjAuNyksIGxhYmVsID0gKCJhLiIpLCBzaXplID04LCBmb250ZmFjZSA9ICJib2xkIikgKwogIGdlb21fdGV4dChhZXMoeCA9MjAsIHkgPSAyMC43KSwgbGFiZWwgPSAoImIuIiksIHNpemUgPTgsIGZvbnRmYWNlID0gImJvbGQiKQoKZ2dzYXZlKEJDX0JBTEFOQ0VEX0dBTV9MTUVSX21lcmdlX2xlZ2VuZF9ncmV5c2NhbGUsIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lID0gIkJDX0JBTEFOQ0VEX0dBTV9MTUVSX21lcmdlX2xlZ2VuZF9ncmV5c2NhbGUucG5nIiwgaGVpZ2h0ID0gOCwgd2lkdGggPSAxNCwgdW5pdHMgPSAiaW4iKQoKI0FMVCBDT0xPUiBCWSBUUkVORApCQ19iYWxhbmNlZF9HQU1fTE1FUl9tZXJnZV9sZWdlbmRfY29sb3JieXRyZW5kIDwtIGdnZHJhdyh4bGltID0gYygwLCA0MC41KSwgeWxpbSA9IGMoMCwgMjEpKSArCiAgICBkcmF3X3Bsb3QocG9pbnRzX3dhdnlsaW5lc19icmF5X2JhbGFuY2VkX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAxLCB5ID0gMSwgd2lkdGggPSAyMCwgaGVpZ2h0ID0gMjApICsKICAgIGRyYXdfcGxvdChCQ19iYWxhbmNlZF9EaXNzaW1pbGFyaXR5X0NvZWZfZXJyb3JiYXJfcmVkdWNlZF9jb2xvcmJ5dHJlbmQgKwogICAgICAgIHRoZW1lKGxlZ2VuZC5rZXkuc2l6ZSA9IHVuaXQoMC41LCAnY20nKSwgI2NoYW5nZSBsZWdlbmQga2V5IHNpemUKICAgICAgICMgbGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTYpLCAjY2hhbmdlIGxlZ2VuZCB0aXRsZSBmb250IHNpemUKICAgICAgICMgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZT0xNCkKICAgICAgICksICNjaGFuZ2UgbGVnZW5kIHRleHQgZm9udCBzaXplKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gMjAsIHkgPSAxLCB3aWR0aCA9IDE5LCBoZWlnaHQgPSAyMCkgKwogICAgZHJhd19wbG90KGdldF9sZWdlbmQoZGlyZWN0aW9uYWxfY2hhbmdlX2xlZ2VuZF9wbG90X2NvbG9yYnl0cmVuZCArIAogICAgICB0aGVtZShsZWdlbmQua2V5LnNpemUgPSB1bml0KDAuNSwgJ2NtJyksICNjaGFuZ2UgbGVnZW5kIGtleSBzaXplCiAgICAgICAgbGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTYpLCAjY2hhbmdlIGxlZ2VuZCB0aXRsZSBmb250IHNpemUKICAgICAgICBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplPTE1KSkpLCAjY2hhbmdlIGxlZ2VuZCB0ZXh0IGZvbnQgc2l6ZSkKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gMjUuNiwgeSA9IDgsIHdpZHRoID0gMywgaGVpZ2h0ID0gMikgKwogIGdlb21fdGV4dChhZXMoeCA9IDIsIHkgPSAyMC43KSwgbGFiZWwgPSAoImEuIiksIHNpemUgPTgsIGZvbnRmYWNlID0gImJvbGQiKSArCiAgZ2VvbV90ZXh0KGFlcyh4ID0yMCwgeSA9IDIwLjcpLCBsYWJlbCA9ICgiYi4iKSwgc2l6ZSA9OCwgZm9udGZhY2UgPSAiYm9sZCIpCgpnZ3NhdmUoQkNfYmFsYW5jZWRfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kX2NvbG9yYnl0cmVuZCwgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSAiQkNfYmFsYW5jZWRfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kX2NvbG9yYnl0cmVuZC5wbmciLCBoZWlnaHQgPSA4LCB3aWR0aCA9IDE0LCB1bml0cyA9ICJpbiIpCgpgYGAKCg==